

%% specify files

matfiles = {'run1a_baseline_sr.mat','run1b_baseline_no_sr.mat',...
    'run5a_crra_sr.mat','run7a_Nstate_sr_fit_v3.mat','run6a_time_varying_s_sr.mat'};

for i = 1:5
    load(matfiles{i})
    stats_list{i} = stats;
end

%% predictability regressions across models, using labor market variables

fprintf('\n\n')

for i = 1:numel(matfiles)
    load(matfiles{i})
    % theta
    fprintf('\n&')
    fprintf('&%2.2f', stats.betas_hpxr_theta_nom)
    % f
    fprintf('&')
    fprintf('&%2.2f', stats.betas_hpxr_fprob_nom)
end